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ABSTRACT 



Context. In the scope of the star formation process, it is unclear how the environment shapes the initial mass function (IMF). While 
observations of open clusters propose a universal picture for the IMF from the substellar domain up to a few solar masses, the young 
association r\ Chamaeleontis presents an apparent lack of low mass objects (m < 0. 1 M©). Another unusual feature of this cluster is 
the absence of wide binaries with a separation > 50 AU. 

Aims. We aim to test whether dynamical evolution alone can reproduce the peculiar properties of the association under the assumption 
of a universal IMF. 

Methods. We use a pure N-body code to simulate the dynamical evolution of the cluster for 10 Myr, and compare the results with 
observations. A wide range of values for the initial parameters are tested (number of systems, typical radius of the density distribution 
and virial ratio) in order to identify the initial state that would most likely lead to observations. In this context we also investigate the 
influence of the initial binary population on the dynamics and the possibility of having a discontinuous single IMF near the transition 
to the brown dwarf regime. We consider as an extreme case an IMF with no low mass systems (m < 0.1 M®). 

Results. The initial configurations cover a wide range of initial density, from 10 2 to 10 8 stars/pc 3 , in virialized, hot and cold dynamical 
state. We do not find any initial state that would evolve from a universal single IMF to fit the observations. Only when starting with a 
truncated IMF without any very low mass systems and no wide binaries, can we reproduce the cluster core properties with a success 
rate of 10% at best. 

Conclusions. Pure dynamical evolution alone cannot explain the observed properties of r\ Chamaeleontis from universal initial con- 
ditions. The lack of brown dwarfs and very low mass stars, and the peculiar binary properties (low binary fraction and lack of wide 
binaries), are probably the result of the star formation process in this association. 



Key words, binaries: general - Stars: luminosity function, mass function 
Chamaeleontis - Stars: kinematics and dynamics - Methods: numerical 
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1. Introduction 

As an imprint of the star formation process and governing the 
evolution of star populations, the initial mass function (IMF) has 
been studied in depth in the solar neighbourhood as well as in 
young open clusters. Particular interest has been devoted to the 
question of the universality of the IMF: is there a unique mass 
distribution resulting from the interplay of physical processes 
of star formation, or does it vary with gas density, metallicity 
dMarks et al.ll2012l) or turbulen ce ? 

Introduced by S alo e ted (Il955l) . the IMF was first described for 
stars in the mass range 0.4 M to 10 M e as a power law 
^(log m) = - °< m~~ r , with T = 1.35 in logarithmic scale. 
Th is field star power-law index was independently established 
bvlKroupa et all {l993a) for 0.5 M to 1 M and extended by 
Massevl J2003l)to 10 Mr; 



F ocussing on close open clu sters, e.g. 
the Pleiades dMoraux et aT|l2003t lLodieu et alj|2007l) IC 4665 
dde Wit et al.ll2006l). a Per dBarrado v Navascues et al.ll2002l) . or 
Blanco 1 dMoraux et al J l2007al) . it was possible to explore the 
system (i.e. incorporating both single objects and unresolved bi- 
naries) mass function in the lower mass regime down to ~ 0.03 
M . Investigations on the shape of the mass function in vari- 
ous environments show some deviations that can be explained 
by uncertainties due to e.g. different sampling, dynamical evo- 



lution, and stellar evolut ion models, but show no evide nce for 
any significant variation (Scalo 2005; Bastian et al. 2010). These 
studies lead to a unive rsal picture of the system IMF down to 
0.03 M (see review of Kroupa et al. 201 1) as a non-monotonic 
function showing a maximum around 0.25 M and a power-law 
tail at the high mass end. Many function al forms can be tailor ed 
to this IMF, e.g. segmented power-laws (Kroupa et al.lll993bl) . a 
log-normal functio n plus a power-law t ail (Chabrier 2003) , or a 
tapered power law (de March iet al.l2005tlParravano et al.B 201 1; 
lMaschbergetl2012l) . 

In this paper, the universality of the IMF is investigated by 
focussing on the dynamical evolution of the stellar g r oup n 
Chamaeleontis. Since its discovery bv iMamaiek et al. (1999) 
this cluster has been the target of many observational studies 
(e.g |Luhmanl2004tlBrandeker et al. 120061: iLvo et al.l2003b. It is a 
young (6-9 Mvr.lLawson & Feigelson 2001; Jilinski et al. 2005), 
close (d ~ 94 pc) and compact group of 18 systems (contained 
in a radius of 0.5 pc). Its system mass function was found to be 
consistent with that of other young open clusters and the field 
dLvo et alj|2004l) in the mass range 0.15-3.8 M , but with a lack 
of lower mass members. This challenges the universal picture of 
the IMF, unless the observed present day mass function has al- 
ready been affected by dynamical evolution. Despite deep and 
wide-field surveys dLuhm an 2004; ISong et al.ll2004t iLvo et al. 
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l2006h . no very low mass systems (m < 0.1 Mp Pb was found 
within 2.6 pc from the center. A recent study by Murphv et al. 
d2010t) reported the discovery of four probable and three possi- 
ble low mass members (m < 0.3 M ) in the outer region, be- 
tween 2.6 and 10 pc from the cluster center. This suggests that 
the lower mass members might have escaped from the cluster 
core due to dynamical encounters and lie at larger radii than 
the more massive members. Moreover, the cluster appears to be 
mass segregated with all the massive stars (m > 1.5 M ) con- 
centrated in its very central region, which supports the picture of 
dynamical evolution. Among the 18 systems 5 are confirmed bi- 
naries and 3 are possible binaries yieldi ng a binary fraction in th e 
range [28%,44%]. As summarised bv iBrandeker et al.l (120061) . 
none of these binaries have a projected separation greater than 
20 AU, and the probability for a star to have a companion at sep- 
arations larger than 30 AU was estimated to be less than 18%. 
This is opposed to t he 58% wide binary pr obability in the TW 
Hydrae association (Brande ker et al.l 12003;). despite its similar 
age. This deficit of wide binaries in rj Chamaeleontis may also 
be explained by their disrupti on through dyn amical interactions. 
In a previous study dMoraux et al.l l2007bl) . we considered 
whether dynamical interactions could explain the lack of very 
low mass systems (m < 0.1 M ) in the cluster core, starting 
with a universal IMF. We applied an inverse time integration 
method by sweeping the parameter space for the initial state in 
order to find those that best lead, as a result of a pure N-body 
simulation, to the observed properties of rj Ch a. This method 
has been applied in nume r ous earlier studies (e.g Kroupal 1 995a: 
IKroupa & Bouvieil 120031: iMarks & Kroupal 12013 to obtain a 
comprehensive picture of the early dynamical evolution of star 
clusters. In our case this was designed as a test of the universal- 
ity of the IMF. A ssuming a log-normal shape for the system IMF 
(Chabrier 2003) we span a large range of initial densities. We 
found that it was possible to reproduce the observations starting 
from a very dense configuration (10 8 stars/pc 3 ) with a success 
rate of 5%. The simulations, however, did not include any pri- 
mordial binaries nor considered the creation of binaries in the 
detailed analysis. The gas was removed initially and we assumed 
that the cluster was in virial equilibrium. 

In the present study, we follow the same method in an at- 
tempt to reproduce the observed state of rj Cha, but we now 
take into account an initial binary population and its evolu- 
tion. In a first set of models, we assume a universal log-normal 
IMF dChabrierl2005l) be fore considering a possible discontinuity 
dThie s & Kroupal l2007l) around the substellar limit, and a trun- 
cated IMF with no system below 0.1 M . The simulations still 
start after the gas has been expelled but virial equilibrium is not 
required. 

The outline of this paper is as follows: we first discuss the 
statistical significance of the deficit of very low mass systems 
(m < 0.1 M ) in the cluster core (Section [2) before describing 
the numerical scheme adopted for the simulations, especially the 
initial conditions and the parameter grid (Section|3). The analy- 
sis procedure is introduced in Section |4] Section |5]presents the 
results obtained when starting with a log-normal IMF. We dis- 
cuss alternative initial conditions in Section |6]before presenting 
our conclusions. 
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1 The lowest mass members have an estimated spectral type around 
M5, leading to masses betw een 0.08 and . 16 M w , depending on the 
adopted evolutionary tracks (Lvo et al. 2003; Luhman & Steeshs 2004) 



Fig. 1. KS test evaluating the similarity of the rj Cha mass distri- 
bution with the log-normal system IMF. The comparison is done 
considering three mass ranges, in order to assess which parts of 
the distribution are alike. We progressively added ten stars in rj 
Cha data set, selected randomly from the IMF within the mass 
range [0.01 ;0. 15] M . Each time a star is added we compute 
a probability from the KS test. This process is repeated for the 
mass ranges [0.15;1] M and [1;4] M . 



2. Statistical issues 

With less than 20 systems in the cluster core, the statistical anal- 
ysis of rj Cha has to be done carefully, especially when consid- 
ering a standard distribution such as the IMF. In the range from 
0.15 to 4 M , the rj Cha mass function (MF) was found to be 
consistent with t he IMF derived f or young embedded clusters 
(iLvo et al. 2004; M ever et al. 2000) and field stars by compar- 
ing the ratio of sta rs with m a ss m > 1 M to stars with mass 
0.1 < m < 1 M . ILvo et al.l d2004l) predicted about 20 mem- 
bers with m < 0. 15 M by comparison with the Trapezium MF. 
None has been found within a 2.6 pc radius despite deep and 
wide searches, indicating a strong deficit of very low mass sys- 
tems in T? Cha. However , using a Kolmogorov-Smirnov (KS) 
test, iLuhman etall d2009l) derived a probability of ~10% that 
rj Cha is drawn from the same IMF as Chameleon I or IC348, 
not revealing significant differences between those distributions. 
One can therefore wonder whether the lack of very low mass 
objects (single stars/brown dwarfs or unresolved binaries with 
m < 0.1 M , hereafter VLMOs) inside a 2.6 pc radius from the 
cluster center represents a significant deviation from the univer- 
sal MF of open clusters. If we ch oose the log-norma l fit to the 
Pleiades system MF as a reference (Moraux et al. 2003), then the 
KS probability for testing the hypothesis that the stellar masses 
in rj Cha are chosen from this log-normal MF is 2.8 %. To assess 
the sensitivity of this result to the data set, we present in Fig. Q] 
the evolution of the KS probability while systems are randomly 
added to the list of known rj Cha members from three differ- 
ent mass ranges ([1-4] M , [0.15-1] M and [0.01-0.15] M ). 
The probability increases uniformly to 80% until eight systems 
with m < 0. 15 M are added, which points out the importance of 
the deficit of VLMOs relatively to medium and high mass stars. 
However the KS probability is already greater than 5% when 
only one such system is added, and we cannot reject the possi- 
bility that this data set is drawn from the MF used as reference. 
As a result, a paucity of VLMOs may be present but it might not 
be statistically inconsistent with the Pleiades MF. 
Nevertheless, even if the deficit of very low mass systems is not 



2 



Christophe Becker et al.: Reverse dynamical evolution of Eta Chamaeleontis 



really significant, it is additional to the other peculiar properties 
that also need to be understood: the lack of wide binaries with 
separation larger than 20 AU, and the presence of mass segrega- 
tion. 



3. Numerical set-up 

In this Section we describe the physical properties that we 
tested in our models: the single initial mass function (of all 
stars counted individually) and the primordial binary properties 
(binary fraction, separation and mass ratio distribution). We 
then review the assumptions corresponding to the lack of gas 
treatment and the density profile, and we present the parameter 
grid that we used for each model: the number of systems, the 
characteristic spatial scale for the density distribution function, 
and the global virial ratio. 



3.1. IMF 

As our main hypothe sis we choose a single IMF of log-normal 
form as sugg ested bv lMiller & Scaiol (1 19791) 



Table 1. Physical properties corresponding to each model. 
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This function was fitted in the 0.1-1 M e mass range to the nearby 
galactic disk MF based on a volume li mited sample w ithin 8 pc 
and yields pL = 0.2 M an d a = 0.55 dChabrieril2005l) . A sim- 
ilar result was obtained by iBochanski et al.1 (1201 Oh for the field 
M-dwarf MF based on a much larger (several million stars), but 
unresolved sample. In order to test its universality, we use the 
single IMF proposed by Chabrier (2005) in our models A, B and 
C. Masses were chosen within a mass range consistent with ob- 
servation, from 0.01 to 4 M©. 

As an alternative, we test the possibility that the single star 
IMF may be discontinuous (model D) with the majority of 
brown dwarfs followi ng their own IMF, as suggested by 
iThies & Kroupal (120071) . Since this may result in a lower num- 
ber of VLMOs in the cluster initially, we might expect this initial 
condition to be more favourable in reproducing the observations. 
We also study the extreme case where the IMF is not universal 
but truncated in the low mass domain, with no system below 0. 1 
M to follow the observations (models E and F). 

3.2. Binary properties 

Observations of young clusters reveal a broad range of binary 

Nu 

fractions bf, defined as bt = — , where Ni, is the number 

f f Nb +N S 

of binaries and N s the number of single objec ts. According to N- 

body numerical simulations (Kroupa 1 995al) . this is consistent 

with a universal initial binary fraction of 100% that decreases 

with time depending on the star cluster density. We thus fixed 

the initial binary fraction to bf =100% for most of our models 

(except model C and F). 

Concerning the mass-ratio distribution, we draw two masses ran- 
domly (in Model A) from the same single IMF in order to get 
the primary and th e secondary masses ; the approach adopted 
by IKroupa d!995al) . For all other mo dels (except Model D) w e 
adopted a flat mass-ratio distribution (Reggiani & Mever 201 ll) . 
and the IMF used to obtain the mass for the primary had to ad- 
justed (see section RTTT i so that once the p airing is done we re- 
trieve the single IMF from Chabrierl (120051) presented above. 





log-normal IMF 
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truncated IMF 
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As for the separatio n distribution, we took the one derived by 
IKroupa et alj d201 lb (see their equ ation 46) to fit the observa- 
tional data for F, G and K field stars (Duquennov & M avorll991l: 
Raghavan et al . 2010) using inverse dynamical population syn- 
thesis and taking into account early evolution of the orbital pa- 
rameters for small period systems. The separation distribution 
ranges from 0.1 AU to 10 5 AU for the models A, B, D and 
E. In accordance with conclusions of the work by Kroupa et al. 
d201l1) . we also chose an initial eccentricity distribution already 
in statistical equilibrium to ensure a thermal final distribution. 
Although observations of multiple systems rather supp ort an es- 
sentially flat distribution of eccentricity (e.g Abt 2006), we note 
that dynamical simulatio ns are hardly se nsitive to the selected 
eccentricity distribution (lKroupalll995bb . We therefore favour 
consistency with previous work to enable direct comparisons. 
Given the observed lack of wide binaries, we consider the pos- 
sibility that no wide binary has formed initially. Instead of as- 
suming a Kroupa-like separation distribution we also explore a 
truncated distribution at large separation together with a smaller 
binary fraction (models C and F, see section l6\2l . 
A summary of all considered models is given in Table Q] 



3.3. Gas 

For all our models, we assume that the gas is already removed 
at the start of the simulations, but that the cluster is not neces- 
sarily relaxed to virial equilibrium. We estimate that our initial 
state depicts a cluster that is between 0. 1 Myr to 3 My r old, de- 
pending on the picture for gas removal dTutukovlll978l) . Indeed, 
there is no clear consensus about the time scale for gas disper- 
sal, estimated from 0. 1 to a few crossing times, depending on the 
mechanism in play (OB star wind, supernovae remnant or stellar 
outflows). In the case of 77 Cha, the mechanism for gas removal 
may involve an ex ternal factor on a ti me scale that could be as 
long as a few Myr. Orte ga et"al] (12009) proposed a common for- 
mation scenario for young clusters in the Scorpio-Centaurus OB 
association (77 Cha, e Cha and Upper Sco), by backtracking bulk 
motions. In this dynamical picture, 77 Cha was bom in a medium 
likely being progressively blown out by strong stellar winds 
coming out from the Lower and Upper Centaums Crux complex. 
Another possibility to expel the gas from t he cluster may involve 
feedback from a massive stellar member. iMoraux et al.l (l2007bl) 
showed that tj Cha initial state might have been very compact, 
with a crossing time of about 2 10 4 yrs (for a total mass of 15 
M Q and a radius of 0.005 pc). In this extreme case, the presence 
of one B8 star (after which the cluster is named) might be suf- 
ficient to remove the gas within 10 5 years. With a cluster age 
estimate of 8 ± 1 Myr, we run the simulations for 10 Myr. 
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Table 2. Parameter grid for model A. We considered all combi- 
nations of N sys , ^pi, and Qi 



Table 3. Selection range adopted for the observational criteria 





20 


30 


40 


50 


60 


70 


Rpi (pc) 


0.3 


0.1 


0.05 


0.03 


0.01 


0.005 


Qi 


0.3 


0.4 


0.5 


0.6 


0.7 





3.4. Density and velocity distribution 

For all models, the systems are distributed spatially using a 
Plummer model 



3M 



sys 



47^ 



[l + (r/Rn) 



21-5/2 



where N sys is the initial number of systems. 
The velocities of each individual object are computed accord- 
ing to this density distribution and to the initial virial ratio 
Qi = E /E pot where E^ is the total kinetic energy of the 
cluster and E po , the gravitational energy. For each model, we 
are thus left with three free parameters: the initial number of 
systems N sys , the Plummer radius Rp\ and the virial ratio, Qj. 



3.5. Parameter grid 

From the shape of the IMF, we estimated an initial value of 
N S ys = 50 by the requirement to have four stars with mass greater 
than 1 M . To cover a wide range of densities at fixed radius, 
we tested N sys from 20 to 70. The initial cluster radius was first 
estimated to fit a constant surfac e density derived from obser- 
vations of star forming regions ( Adam s"et al.l 2006|). giving 0.3 
to 1.0 pc for 50 systems. The study of iMoraux et al.l (l2007bl) 
showed that a dense initial configuration was necessary in order 
to eject enough members from the cluster core and reproduce 
the lack of VLMOs. To favour dynamical interactions, we took 
a radius varying from 0.3 to 0.005 pc, yielding a density range 
from 500 stars/pc 3 to 10 8 stars/pc 3 . In order to assess the effect 
of initial equilibrium, we tested cold, gravitation-dominated con- 
figuration (Qi = 0.3) and hot, initially expanding configuration 
(Qi = 0.7). 

For each model (A to F), an initial configuration is characterized 
by a combination of {N sys , Rp\, Qi} from the values given in 
Table[2] In total, 180 arrangements were tested for each model. 



Criterion 


Range 


Restricted to 


o V SLCIIla 


a/, _ [14. 971 


/ O.J pc 


Massive stars 


N 2 = [2A] 


1 \J.J pc 

m > 1.5 M 


VLMOs 


N 3 = [0, 1] 


1 \ jL.KJ pc 
m < 0. 1 M:) 


Halo 


N 4 = [0, 1] 


0.5 < r< 10 pc 
m > 0.5 M Q 


Binary fraction 


N 5 = [22,50]% 


r < 0.5 pc 


Wide binaries 


N 6 = [0, 1] 


[50; 400] AU 
r < 0.5 pc 


Time 


[5,8] Myr 





defines its tidal radius: 



GM 



4A(a-b; 



1/3 



1.4 M 1/3 pc 



where M is the cluster t otal mass (i n solar masses) and A and B 
are the Oort constants (Kind fT962l) . Given the parameter grid, 
the initial r t value varies between 3.1 and 4.8 pc (the estimated 
valu e for 7] Cha is around 3.5 pc assuming a total mass of 15 
M^. lLvo et 311(12004 1. Objects are considered as being ejected 
and then removed from the simulation as soon as they are 
further than twice the cluster tidal radius from the center. 
For each initial configuration {N syi , R?\, Qi} we generated 
200 simulations, changing only the random seed, for statistical 
purposes. Every simulation computed the cluster dynamical 
evolution for 10 Myr. 



4. Analysis procedure 

In order to retrieve as much information as possible we anal- 
yse our set of simulations in two different ways f or each model. 
First we consider the same analysis procedure as in M oraux et al.l 
(l2007bl) that aims at finding final states that fit the observational 
data. Secondly, in order to better understand the results of the 
first analysis, we perform a statistical analysis. Both methods are 
based upon a set of constraints derived from the observations. 



3.6. N-body code 

We use the NBODY3 code (lAarsethl 19991) that performs a direct 
force summation to compute the dynamical evolution of the clus- 
ter. Close encounters are treated b y Kustaanheimo-Stiefel (KS ) 
regularization for hard binaries (iKustaanheimo & Stie fel 1965), 
which uses a space-time transformation to remove the singular- 
ity and then si mplify the two body treatm ent, or chain regular- 
ization method (Mikkola & Aarseth 1990) for few body interac- 
tions (e.g. binary-single star). There is no stellar evolution. 



4.1. Observational criteria 

We use a set of criteria described below to evaluate if a simu- 
lation at a given time is close to reproducing the observations. 
Each criterion is associated with a range of validity assuming 
Poisson statistics: a criterion i is satisfied if A^,- € [<9; — \/0~i, Oi + 
\/Oi], where Af, is obtained by simulation and (9, is given by 
the observations. A summary of the chosen ranges is given in 
Tabled 



3. 7. Modelling procedure 

The time evolution of the initial conditions described earlier pro- 
duces output of positions and velocities for each star every 0.05 
Myr. The NBODY3 output files also provide details for close 
binaries (semi-major axis, eccentricity) identified as bound dou- 
ble systems. The stellar cluster is put in a galactic potential that 



Number of systems (A^i) To account for the membership and 
compactness of the core, we consider the total number of sys- 
tems in a 0.5 pc sphere. Since 18 systems have been observed 
within the core radius, we choose the range of 14 to 22 systems 
for a simulation to fulfil this criterion. Unless mentioned other- 
wise the term system refers to a single object or a binary of any 
mass within the stellar or substellar domain. 
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To take into account observational limitations in the compar- 
ison between simulations and observations, we identify bina- 
ries as closest neighbour pairs in projection (i.e. not neces- 
sarily bound) with a separation smaller than 400 AU (which 
corresponds to 4" at the cluster's distance). At larger separa- 
tions binaries are observationa lly identified as two single objects 
dKohler & Petr-Gotzensll2002l) . 

Number of massive stars (A^) Since three systems were found 
in the central region with a mass m > 1.5 M & , we require to 
have between 2 and 4 of them in the simulations. When counting 
massive stars, a binary system is considered as a single object 
with a mass corresponding to its total mass. 

Number of systems in the halo (A^) No potential cluster mem- 
ber has been identified by the ROSAT All Sky Survey (sensitive 
to late-K type stars) outside the cluster core up to a distance of 
10 pc. This translates into the following criterion: less than one 
cluster member more massive than 0.5 M E , must lie within the 
distance range [0.5-101 p c from the cluster center. 
Recently Mur phy et al.l (1201 Oh have discovered four probable 
and three possible less massive members (in the spectral range 
K7 to M4, i.e. O.1M < m < 0.3 M ) at a distance between 2.6 
pc and 10 pc from the cluster centre. However, since the status 
of these candidates is not confirmed, we will check a posteriori 
that some simulations matching all other criteria do produce a 
number of low-mass halo stars that is consistent with the small 
number suggested by Murphy's study. 

Number of very low mass objects (A4) No system with m < 
0.1 M q has been fo und within 2.6 pc radius from the cluster 
centre (Luhman 2004). The associated criterion is to have either 
zero or one of this kind of object left in the simulation. 
The absence of very low-mass systems is observed for both sin- 
gle objects and companions at a separation larger than 50 AU. 
In our simulations, the number of VLMOs is therefore the total 
number of companions (within a separation range of [50-400] 
AU), single objects and close binaries (separation smaller than 
50 AU) whose mass is below 0.1 M Q . 

Binary fraction (JV5) iBrandeker et al.l (120061) identified 5 bina- 
ries and 3 candidates for a total of 18 systems in the core re- 
gion. Considering the average value of 6.5 binaries, this gives 
an observed binary fraction of 36% and the validity range for 
this criterion is from 22% to 50%. Since binaries wider than 400 
AU are considered as two separate single stars in our analysis, 
the simulated binary fraction is already of the order of 50% be- 
fore any dynamical evolution for models A, B, D and E because 
of the initial period distribution. Therefore this criterion is not 
expected to be critical. 

Number of wide binaries (Afc) f] Cha does not contain any bi- 
nary with a projected separation greater than 30 AU. This was 
put into the following constraint : we require the model not to 
contain any binary with a separation larger than 50 AU. We 
choose a loose cut on separation to be more conservative and 
to take projection effects into account. In the following we refer 
to the number of wide binaries as the number of binaries with 
separations larger than 50 AU and smaller than 400 AU. 



Age With an initial state estimated to be between 0.1 to 3 Myr 
(see section \33[ , and an age for 77 Cha taken to be 8 Myr, we 
require the simulations to be in the age range from 5 to 8 Myr. 
We also require the time window during which the other criteria 
are fulfilled not to be smaller than 1 Myr, to exclude transient 
states. 

4.2. Probability maps 

Since it appears very difficult to satisfy all criteria simultane- 
ously for most models, we refine our analysis and build a prob- 
ability to estimate how likely a set of simulations reproduces 
each observational constraint independently. At each time step 
we compute the probability «,(?) for the simulation to fulfil a 
criterion i. This probability is calculated from the normalized 
histogram generated from the simulations by summing all the 
bins in the range [Oi — \fOi 1 Oj + \fO^\ ■ Statistical scatter is dealt 
with using a smoothed histogram in case of a poor bin sampling. 
If none of the 200 simulations recover the range associated to 
the observed value, we set the probability to 1/200, regardless of 
the gap separating this interval to the value of the first non-zero 
bin. In the case of a complete mismatch between observation and 
model, this method does not provide more information than an 
upper limit. 

The probability a,(f) can be calculated for each configuration 
{A^sys, Rph Qi} an d eacn model. In particular we can produce 
maps of CLi(ti jtn ) in coordinates of Af sys and Rp\ for a given Qi 
and a given model (see e.g. Fig. |4j, where f,-. m corresponds to the 
time, in the range [5,8] Myr, at which a, (f ) is maximum. 

5. Results from our standard model (model A) 

In this section, we discuss the results given by model A to test 
whether rj Cha can be reproduced from a universal log-normal 
single star IMF, with 100% binary fraction and random pairing. 
This model is a first guess, based on standard assumptions. The 
analysis presented below motivated us to relax some assump- 
tions (section|6]l. 

5.1. Reproducing rj Cha 

The criteria described in the previous section allow us, when 
used together, to check the ability of model A to reproduce the 
observations for a given set of initial parameters. Considering 
each of the 200 realizations for all configurations {N sys , Rp\, Qi}, 
we apply these criteria at each time snapshot to see if they can 
all be satisfied simultaneously. Table. [4] shows a summary of this 
procedure for a specific value of virial ratio (Qi = 0.5) and num- 
ber of systems (N sys = 20) for the first 4 criteria (thus without 
any constraint on the binary properties nor the age). 
Even if most runs satisfy the first criterion on the number of 
systems, this is valid only for a given time range, in which the 
next criterion will have to be fulfilled. The most important result 
is that the percentage of runs passing the selection drops to zero 
as we apply the fourth condition on the number of VLMOs for 
all the initial configurations with N sys > 20, and to 0.5% at best 
forays =20. 

This indicates that this criterion is very difficult to fulfil simul- 
taneously with the other criteria. A dense initial state is nec- 
essary to remove all (or almost all) very low mass members 
from the cluster by enhancing two-body encounters, especially 
as more objects are released during the processing of wide bina- 
ries. However this tends to quickly inflate the inner core, acting 
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Q=0.5, N=70, R=0.005 pc Q=0.5, N=20, R=0.3 pc 




Time (Myr) Time (Myr) 

Fig. 2. The left panel displays the velocity of ejected members at the time of ejection (blue asterisks for single objects, and red 
open circles for binaries) for a dense configuration from model A, with N sys = 70, Qi = 0.5, and Rp\ = 0.005 pc. The right panel 
corresponds to a sparse configuration with N sys = 20, Qi = 0.5, and R?\ — 0.3 pc). The half mass radius evolution is superimposed 
(thick line), as a footprint of the stellar density for both cases. 



Table 4. Results of the quantitative analysis for different values 
of Rp\, and for N sys = 20 and Qi = 0.5 for model A. For each 
value of R p i we apply the four first criteria one after the other. 
Every time a criterion is added, we compute the number of runs 
that fulfil the condition. As a result of our successive elimina- 
tion scheme, only three simulations satisfies the first four crite- 
ria. However, none of those fulfil all six criteria simultaneously. 



Criterion 


0.005 


0.01 


0.03 


0.05 


0.1 


0.3 


Systems 


195 


199 


197 


199 


199 


200 


+ Massive stars 


96 


114 


139 


150 


154 


143 


+ Halo 


21 


41 


67 


89 


100 


72 


+ VLMOs 





1 


1 





1 






in opposition to the criteria on the number of systems (A^) and 
massive stars (A/2) and increasing the number of stars in the halo 
(A^). This is illustrated for the dense initial configuration with 
{A^ys = 70, R Pi = 0.005 pc, Qi = 0.5} in the left panel Fig. 
There is a peak of cluster members ejection^ before 1 Myr with 
velocities as high as 60 km/s (especially for single objects re- 
leased by binary decay). As an imprint of this highly dynamic 
phase the cluster undergoes a fast expansion phase, shown by 
the increase of the half mass radius from 0.01 pc to 1 pc within 
1 Myr. Once the density has fallen off, the dynamics involves 
softer interactions (secular evolution) and the number of ejected 
members decreases along with their velocity. During this phase 
the cluster expands slowly until reaching virial equilibrium. 
As a result of the fast relaxation phase the VLMOs are ejected 
efficiently but the numbers of systems (Ni) and massive stars 
(Ay remaining in the cluster core are too small. In addition, the 
core expansion adds many solar-type systems to the halo, incom- 
patible with the criterion A3. We can move to a less dense initial 
state to try to improve the results, but then the expansion is too 
slow and the number of VLMOs inside a 2.6 pc radius (A/4) re- 
mains almost constant with time. When starting with a sparse 
configuration (A^ sys = 20, Rp\ = 0.3 pc and Qj = 0.5 ), we do not 



ejected member are any object unbound to the cluster and being at 
a distance larger than twice the half mass radius from the cluster center 



see any peak of ejection at earlier times, and the half mass ra- 
dius increases slowly and linearly in time (Fig. [2] right panel). It 
seems therefore that a compromise on the initial density has to 
be found in order to eject most of the VLMOs while retaining 
a dense enough core (compatible with criteria N\ and A/2) and 
without populating the halo. 

To better understand the cluster dynamical evolution, we show 
in Fig. [3] the evolution of the six quantities constrained by the 
observations for the 200 realizations that started with an inter- 
mediate density (N sys = 40, Qi = 0.5 and R P \ = 0.05 pc). The 
range corresponding to each criterion is delimited by solid lines 
in each panel. 

First, it is interesting to note from the top left panel that the num- 
ber of systems does not actually start at the setup value 40, but 
around 53 in average. This is mainly due to not counting bound 
pairs with separations larger than 400 AU as binaries but as two 
single objects, thus increasing the total number of systems. This 
can also be seen in the lower middle panel, where the binary 
fraction is initially around 46%, instead of 100% as set up. Then, 
during the cluster early evolution phase, binaries are processed 
more or less efficiently due to dynamical interactions, depend- 
ing on their separation and on the initial density. In this case the 
binary fraction decreases from 46% to 43% within 0.5 Myr. As 
a consequence of the binary disruption the number of systems 
inside the inner core increases slightly during the first 0.5 Myr. 
After this phase, dynamical interactions are softer. The secular 
evolution tends to inflate the core, slowly dispersing the cluster 
members, decreasing the number of systems (N\ ) and wide bi- 
naries (A/5) in the inner core, and increasing the number of stars 
in the halo (A/3). 

The number of VLMOs inside a 2.6 pc radius (A/4) evolves in a 
similar way to the number of systems in the inner core (see bot- 
tom left panel). Note that the number of very low mass systems 
expected from the initial conditions (log-normal single star IMF, 
100% binary fraction and random pairing) should be around 6 
for N sys = 40. However, the number of VLMOs is already ~ 13 
at t = Myr (bottom left panel of Fig. [3j, due to the fact that any 
very low mass (m < 0. 1 M ) companion at separation largerthan 
400 AU is counted as a single object. This number remains con- 
stant during the first 0.5 Myr as binary disruption compensates 
for the ejection process. Then, later in the cluster evolution, the 
number of VLMOs (A/4) decreases slowly, as does the number 
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Fig. 3. 3D histograms showing the evolution of (a) the number of systems in the core N\ , (b) the number of massive stars in the core 
N2, (c) the number of systems in the halo N3, (d) the number of VLMOs inside a 2.6 pc radius N4, (e) the binary fraction in the core 
N5, and (f) the number of wide (separation greater than 50 AU) binaries N$ as a function of time for model A and for the initial 
configuration {N sys = 40, Rp\ = 0.05 pc, Qj = 0.5}. The histogram density corresponds to the number of simulations (out of 200) 
that fall in a given bin. The dashed and solid red lines correspond to the observed value and the acceptable range respectively. 



of systems. However, it remains larger than five for most of the 
simulations starting with N S y S = 40, Rp\ = 0.05 pc and Qi = 0.5. 
Overall, the success rate for the simulations to reproduce the ob- 
servations is zero for all the initial configurations in model A. 

5.2. Best-fitting initial conditions 

In order to know which criterion is the most stringent and how 
the model hypothesis could be modified to reproduce 77 Cha, we 
performed the analysis based on probability maps of «,(?,■.,„) de- 
scribed in section 1431 Fig. |4] reveals the regions (in {7Y sys , ^pi} 
coordinates, for Qi — 0.5) which most likely satisfy a given ob- 
servational criterion. Below we review the inability of the simu- 
lations performed in model A to reproduce the observations, in 
the light of Fig. |4] Note that the time constraint is not mentioned, 
but was applied to produce all the discussed probability maps. 

A^i We notice that the probability to fulfil the criterion on the 
number of systems in the core drops if we start with large initial 
values for N sys and Rp\ since the density becomes too small 
to remove enough systems from the inner core by dynamical 
interactions. On the other hand, when starting with a low N sys 
and Rp\, the number of systems that remain in the core becomes 
rapidly too small. 



N2 The criterion on the number of massive stars seems to be 
easy to reproduce and does not strongly depend on the initial 
parameters although there is a small trend in favour of less 
dense cases or large value of N sys . 



N3 Considering the number of systems in the halo, it is clear 
that this criterion is best matched with the smallest N sys because 
less objects can be ejected in the halo. For N sys = 30 or 40, this 
criterion is more easily fulfilled for either large Rp\ (as lower 
density leads to fewer ejections), or small Rp] (as high density 
induces fast ejections, leading to large projected distances by 5 
Myr). Intermediate values of Rp\ result in too many slow-moving 
ejected stars that will remain in the vicinity of the cluster. For 
N sys > 40 the criterion on (A^) is very badly reproduced for any 
value of Rp\. 



N4 The result for the VLMOs is important since the region of 
agreement is very narrow. This shows that this criterion together 
with N3, is the most stringent. It requires a low value of N sys to 
minimize the initial number of VLMOs to eject, and a low Rp\ 
to maximize the dynamical encounters and eject these objects 
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Fig. 4. Summary for all tested configurations for model A with a virial ratio Qi = 0.5. For convenience the value for the quality 
measure, a,(f,, m ) is also indicated by the dot size, with a logarithmic scale. For each criterion, the colour (and size) of the circles 
indicates the probability for the simulations to reproduce the observations. The results obtained with a different virial ratio Qi are 
very similar. 



efficiently. 



N5 and Ns We notice that the map for the binary fraction (N$) 
does not indicate a large dependence on the parameters with an 
overall good agreement with the observations. The separation 
map (A^) reveals a higher probability for the dense cases, which 
process the widest binaries and expand fast enough so that less 
binaries are present in the central region. 



Although very narrow, the overlap region between the agree- 
ment maps of the various criteria (especially those for N3 and 
N4) seems to indicate that suitable initial configurations may be 
found (e.g see Fig. |4]l for intermediate to low N sys and low Rp\ 
(except for the lowest values for which the criterion on the num- 
ber of systems N\ is not well fulfilled). However this result is 
misleading for the criteria are not independent. There is a sig- 
nificant anti-correlation between the number of stars in a 10 pc 
halo (N3) and the number of VLMOs in a 2.6 pc radius (N4) at 
low and medium initial densities (as seen Fig. [5] for N sys = 30 
and Rpi = 0.05 pc). In these cases the constraints on N3 and N4 
tend not to be compatible. At higher densities and especially for 
Rpi = 0.005 pc, this anti-correlation becomes negligible. Both N3 
and N4 get very small: the strong dynamical interactions remove 



all VLMOs from the core, and most ejected objects travel much 
further away than 10 pc within 5 Myr due to the high ejection ve- 
locities. However, the dynamical interactions are so strong that 
it is very difficult to retain anything in the cluster core and the 
number of systems N\ becomes too small. Fig.[6]shows the cor- 
relation between N\ and N4 at t = 7 Myr for the 200 simulations 
starting with N sys = 30 and Rp\ = 0.05 pc. In all cases, when the 
number of VLMOs is less or equal to one, the number of sys- 
tems in the core is smaller than 10, making these two criteria 
incompatible. We found a similar correlation for all the initial 
configurations tested by our simulations. The negative result of 
the previous analysis indicates that ejecting all VLMOs from the 
cluster core and keeping enough systems in a 0.5 pc sphere is a 
major challenge. 

The results obtained for different values of Q\ (from 0.3 to 0.7) 
cannot be statistically distinguished from those obtained for the 
state initially in virial equilibrium (Qi = 0.5). 

5.3. Summary and comparison with the previous study 

Using standard initial conditions corresponding to model A, we 
tested many different configurations, varying the density and 
virial ratio. The main conclusions from this analysis are that 

- starting with a single log-normal IMF with a peak value 
jj. = 0.2 Mj, and a deviation a = 0.55, and assuming an ini- 
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Fig. 5. Correlation map between the number of VLMOs (JV3) and 
the number of solar-type stars in the halo (N4) at t = 7 Myr for 
N sys = 30, R Pi = 0.05 pc and Q t = 0.5 (model A). 
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Fig. 6. Correlation map between the number of systems (A^ ) and 
the number of VLMOs (N3) at t = 7 Myr for JV sys = 30, Rm = 
0.05 pc and Q t = 0.5 (model A). 



tial binary fraction of 100%, random pairing and a Kroupa- 
like period distribution for the binary population, does not 
allow the simulations to reproduce the observations for any 
configuration {N sys , R?\, Qt}; 

there is no hint of an improvement at the edge of the param- 
eter grid, suggesting that our failure to find a solution is not 
a consequence of using a limited parameter space. 



In iMoraux et al . (2007b) the best fitting set of initial parame- 
ters gave a success rate of about 5%, whereas in our analy- 
sis of model A it is 0%. The apparent divergence between our 
results and Moraux's is the consequence of the initial condi- 
tions. In the previous study, the chosen IMF corresponded to the 
system mass function obtained after binary processing (as it is 
observed in the field or in clusters), and binaries were consid- 



ered as unbreakable objects, unable to exchange energy to the 
cluster by modifying their orbital properties. Here, the system 
IMF peaks at higher masses, generating more systems with mass 
m > 0.5 M initially, potentially increasing the number of them 
that could end up in the halo. This makes the criterion on N$ 
more difficult to fulfil in the present study. Besides, binary dis- 
ruption can significantly alter our ability to reproduce criterion 
N4. Even though there are less very low mass systems initially, 
many objects with m < 0.1 M belong to a binary system with 
a separation larger than 50 AU or have been be released by bi- 
nary decay. In both cases, these objects will be accounted for in 
the number of VLMOs (A^), and this criterion is therefore not 
improved. 



6. Alternative initial conditions 

We discussed above the importance of the binary population in 
shaping the system IMF and hosting VLMOs that may be re- 
leased in the cluster core. Since these processes depend strongly 
on the binary properties (mass ratio and separation distribu- 
tions), we will now describe how they may be adjusted (model B 
and C) to better reproduce the observations. We will also discuss 
the possibility that the single star IMF might be discontinuous 
around the substellar limit (model D), which may help to reduce 
the initial number of VLMOs in the core. We will then present 
the results obtained in the extreme case when starting with an 
IMF truncated at 0.1 M and a binary fraction of 100% (model 
E) or less (model F). 



6.1. Binary pairing (Model B) 

In model A, we chose for simplicity to pair binaries randomly 
from the same singl e IMF Neverthe l ess, recent studies of both 
the g alactic field (Raghav an et alJ 



2010; Reggiani & Meve 



5 



1201 ll) and star forming regions (iRraus et all 120081 12011 
indicate that, whereas there is no clear and unique best fit, a flat 
mass ratio distribution may be a better fit than a random pairing. 
Since this would result in a slightly smaller number of very 
low mass companions, we may expect the criterion on N4 to 
be better fulfilled. To implement it in our initial conditions for 
model B, we sample the primary mass from a primary IMF and 
then draw the secondary mass according to a flat mass ratio 
distribution. This requires a slight change in the parameters of 
the primary IMF, in order to reproduce the log-normal single 
IMF. This gives pi = 0.32 M and a = 0.55 (instead of 0.2 M 
and 0.55 in case of random pairing). The binary fraction and 
separation distribution are the same as in model A. 
We ran simulations for Qt = 0.5 only, A^ ys =20 and 7V sys =40, 
with the same range for Rp/ as before. Results show that the 
agreement probability on the number of VLMOs (04) is larger 
(by a factor of two to three), contrary to the probability 03 on 
the halo that is smaller (by a factor of about 2). This is due to the 
shift towards higher masses of the primary MF yielding more 
objects with m > 0.5 M that may end up in the halo. Overall, 
no significant improvement is observed when using a flat mass 
ratio distribution since again no simulation is able to reproduce 
the observations. 



6.2. Separation distribution (Model C) 

In the following, we discuss the possibility that no wide binary 
formed initially by assuming an initial separation distribution 
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Fig. 7. Cumulative separation distribution function (solid line) 
obtained from the Kroupa period distribution (see section |3j. 
This highlights how the distribution is truncated in Model C in 
the case of ft, = 0.8. Expression Q] yields = 0.6, which im- 
poses to truncate the separation distribution to 370 AU, so that 
Nhb = 0.6 x Nb- 



Table 5. Binary fraction and separation cut-off for model C. 



fb (%) 


100 


90 


80 


70 


60 


50 


40 


Cut-off (AU) 


730 


570 


370 


240 


150 


90 


50 



similar to model A but truncated at large separation. The cut- 
off separation value and the initial binary fraction are linked to 
each other and we explain below how they can be evaluated pro- 
viding the final binary fraction. 

Following the simplistic argument that all binaries with a separa- 
tion smaller than a given value (hard binaries) survive through- 
out the simulation, and that any wider binary is destroyed, we 
can express the initial binary fraction fb in terms of the initial 
hard binary fraction = Nhb/Nb (where Nb is the number of 
binaries and Nbi, the number of hard binaries) and the final bi- 
nary fraction //: 

/' 

fh ^ {l+fblfkb-fb 

Taking a final binary fraction of 36% given by the observations, 
we consider different values of fb, ranging from 36% to 100%. 
The corresponding initial hard binary fraction ranges from 100% 
(all binaries survive) to 53% (about half the binaries survive) re- 
spectively. To follow the observations, we identify as hard bi- 
naries (that will not be destroyed) those with separations lower 
than 50 AU. From the initial hard binary fraction, we then es- 
timate t he corresponding separation cut-off, assuming a Kroupa 
(Il995bl) distribution below this value. For example, we need a 
cut-off at 730 AU for a hard binary fraction of 53%. The lowest 
possible cut corresponds naturally to 50 AU, to get 100% hard 
binaries. Fig. [7] illustrates this process in the case of an initial 
binary fraction of fi, = 0.8, which gives fbb = 0.6. This initial 
hard binary fraction is obtained when applying a cut-off in the 
separation distribution around 370 AU. 

We can wonder why these binary properties would result from 
the cluster formation process and this needs to be compared 



to what is observed in star forming regions. In dense envi- 
ronment such as the Trapezium fh ~ 60%, whereas in sparse 
regions like Taurus ~ 90% dDuchend [l999: Kirk & Myers 
2012). A plausible explanation for this difference is that all star 
forming regions start their evolution with a high binary frac- 
tion and the wide binaries are fu rther disrupted in dense en- 
vironments within 1 Myr (see e.g iMarks & Krou"pall2012l) . In 
the Trapezium, very few bin aries with a separ ation larger than 
1000 AU have been found dScallv et alJI 1999b which supports 
this idea. For instance it was possible to reproduce the evo lution 
of the ONC dKroupa et al.ll2Q01b IMarks & Kroupall2012l) start- 
ing with 100% binaries and a density of 10 5 stars/pc 3 , includ- 
ing the deficit of [200-500] AU binaries compared to th e sepa- 
ration distribution for field binaries from Raghavan et al. I (l2010l) 
(Reip urth et al.l2007h . In our simulations starting with N sys = 20 
and Rpi = 0.05 pc, the initial density is very similar (2 x 10 5 
stars/pc 3 ). However, the adopted separation cut-off at 50 or 90 
AU cannot be explained by dynamical encounters since these 
separation limits are much lower than the initial mean neighbour 
distance (around 2200 AU). Nevertheless, it is still possible that 
the binary fraction may be set up during the formation process 
and/or during the gas-rich phase which is not covered in our sim- 
ulations. 

We ran the simulations for <2, = 0.5, N sys = 20 to 70, and 
Rpi = 0.005 pc to 0.3 pc (model C). The parameters used for 
the binary fraction and separation cut-off are given in Table [5] 
A fiat mass ratio distribution (as in model B) has been used to 
generate the secondary masses. 

Fig.[8]shows the evolution of the probability in the case 

A^ S y S = 20 and Rpi = 0.1 pc for the six criteria as a function of 
the adopted separation cut-off. For a large separation cut-off, the 
probability of agreement for the binary fraction is low (< 0.2). 
This is worse than what was obtained for models A and B, for 
which no cut-off was applied to the Kroupa-like separation dis- 
tribution. This is because (1) more binaries have a separation 
lower than 400 AU and are thus identified as binaries in the anal- 
ysis procedure leading to a higher initial binary fraction, and (2) 
the high binary fraction remains almost constant in time, unless 
the initial density is very high. An improvement is naturally seen 
for the criteria on the binary fraction as well as on the number 
of wide binaries when the separation limit gets smaller (< 100 
AU). Applying a cut-off at 50 AU corresponds to removing the 
constraints on the binary population since we already start with 
what is observed (no wide binaries and fb = 40%). The proba- 
bility «3 for the halo is also increasing, from 0.08 to 0.4 for the 
lowest cut-off. For the number of systems {N\), the number of 
massive stars (Ni), and the number of VLMOs (Na) the prob- 
ability does not change significantly. This may be surprising at 
first, especially for N4, as less VLMOs will be produced by bi- 
nary decay. However, this effect is compensated by the slower 
dynamics making it more difficult to eject the VLMOs from the 
core even though they are less numerous. 
Nevertheless, the analysis reveals two configurations (Q; = 0.5, 
A^ys = 20, and Rpi = 0.05 pc and 0.1 pc) for which some simu- 
lations satisfy all criteria if the separation cut-off is 50 AU. We 
found respectively one and three simulations out of 200 that ful- 
fil all the observational constraints. 

To check whether t hese successful run s are consistent with the 
recent results from Mur phy et al.l (1201 Ol) . we look at the number 
of low mass systems in the mass range [0.08, 0.3] M located at 
a distance range [2.6, 10] pc from the cluster centre. We find be- 
tween zero to one of these systems, which is possibly too small 
compared to the detection of four probable plus three possible 
candidates. 
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Fig. 8. Evolution of the agreement probability for the different 
criteria, when setting a maximum value for the binary separa- 
tion, varying from 50 to 730 AU, and adopting a fiat mass-ratio 
distribution (model C). The results are given for {(2,-0.5, N sys = 
20, R p[ = 0.1 pc}. 

6.3. Treating brown dwarfs as a separate population (Model 
D) 

So far we have considered a continuous IMF that extends to 
the su bstellar regime (down to 0.01 M ), but iThies & Kr oupa 
(2007) suggest that brown dwarf (BD) formation may be differ- 
ent to star formation (based on their binary properties), which 
would lead to a discontinuous mass distribution for single ob- 
jects. This assertion is still a matter of debate, but nonethe- 
less finds observatio nal supporQ from the mass function of 
young s tar clusters (IThies & Kroupal 12008}) and BD binaries 
surveys dKraus & Hillenbrandl2012l) . lParker & Goodwin! (1201 lb 
excluded pure dynamical evolution as a possible explanation for 
the observed differences between the separation distributions of 
stellar and substellar binaries, implying that it may be a pristine 
feature (or set during the very early evolution). From a theoreti- 
cal point of view, the process of BD formation remains unclear 
and may involve a star-like collapse within a turbulent medium 
(e.g. lWhitworth & Stama tellos 20Qa) or a more specific channel 
of early ejection of gas eous clumps dReipurth & Clarke! 1200 it 
iBasu & Vorobvovll2012l) . O ther plausible mechanis ms suggest 
massive disc fragmentation (Stamatell os et al.ll2007l) or gravita- 
tional instabilities i nduced in disks a s a result of encounters in 
embedded clusters dThies et aill2010l) . 

To evaluate the possibility that the single IMF may be discon- 
tinuous (model D) , we consider initial con ditions corresponding 
to the results from lThies & Kr oupa (2008). We adopted two log- 
normal single mass functions with jj. = 0.08 M@ and o = 0.69, 
but one corresponds to stars and is limited to the mass range 
[0.07, 4] M and the other one corresponds to BDs and very 
low mass stars (VLMS) with 0.01 < m < 0.15 M Q . There is 
an overlap between the two mass functions in order to end up 
with a continuous system IMF consistent with the universal pic- 
ture of the IMF. Each population (stars, and BDs + VLMSs) is 

3 A recent review (Jeffries 2012) emphasize that the lack of coher- 
ence and completeness of the observations do not allow firm conclu- 
sions. 



treated separately and the BDs and VLMSs to stars ratio is as- 
sumed to be 1/5. The BD and VLMS binary fraction and the 
star binary fraction are respect ively 30% and 100% and there are 
no mixed BD/VMLs binaries (Kroupa et alJl201 lb . For simplic- 
ity we generate the binaries for each population using random 
pairing and the same period distribution with no separation cut- 
off. The latter hypothesis is not realistic since fi eld BD binaries 
are kn own to have a tighter period distribution (Burgasse r et alj 
l2007h that cannot be exp lained by pure dynamical evolution 
(Park er & Goodwinll201 ll) . Nevertheless this will have a very 
limited impact on our results, since the number of BD binaries 
is one or two in average (if starting respectively with N sys = 20 
or N sys =40). 

We ran simulations in the virialized case for N sys — 20 and Rp\ 
within {0.05, 0.1} pc as well as for 7V sys = 40 and R Pi = 0.05 
pc. As a result of the analysis the improvement over our pre- 
vious simulations is limited: no simulation matches all obser- 
vations of the r\ Cha cluster. Compared to the standard case 
(model A), the main improvement lies in the probability aj, to 
reproduce the halo, which is mainly explained by the shift to- 
wards lower masses of the system IMF. At best the probabil- 
ity increases from 0.29 to 0.55 in the case with iV S y S = 20 and 
Rpi =0.1 pc. Note that this probability is also higher than when 
we applied a separation cut-off at 50 AU (model C, see Fig|8j. 
However, the probability of agreement for the criterion on the 
VLMOs decreased compared to the result from model A (from 
0.03 to 0.005 for N sys = 20 and R n = 0.1 pc). We can understand 
this by counting the mean number of VLMOs at t=0.6 Myr, af- 
ter the binary breaking phase: for N sys = 20 and Rp\ = 0. 1 pc we 
find Nvlmo = 8, compared to Nvlmo = 6-5 in the standard case 
(model A) and to Nvlmo — 4.2 in the extreme case of model C 
starting with 40% binaries and a cut-off at 50 AU in separation. 
Since the number of VLMOs is a strong constraint, this compar- 
ison shows the limited effect of the changes adopted for the BD 
population. 

6.4. Truncated IMF at the low mass limit and truncated 
separation distribution (Models E and F) 

The previous analysis indicates that the observational result re- 
garding the number of VLMOs in 77 Cha is particularly diffi- 
cult to reconcile with the other constraints, in particular with the 
number of systems in the core and the absence of solar-type star 
in the halo (see Fig. [5]and[6]l. In the following we consider the 
extreme scenario where the IMF is not universal and no very 
low mass system (m < 0.1 M Q ) has formed initially. To do so, 
we generate primary masses from the same primary IMF (peak- 
ing at 0.3 M©) as in model B, but truncated at 0.1 M e , and use 
a fiat mass ratio distribution (without any truncation on the sec- 
ondary mass) 

We first ran simulations starting with a virialized Plummer 
sphere with 100% binaries drawn from the Kroupa separation 
distribution (model E). iV S y S and Rp\ are chosen within {20, 30, 
40, 50, 60} and {0.005, 0.01, 0.03, 0.05, 0.1, 0.3} pc respec- 
tively. We discarded the larger value N sys = 70 since it would 
give a cluster starting with too many systems to fulfil the criteria 
on the number of systems in the inner core without populating 
the halo. As a result no simulation fulfilled all 6 criteria. This as 
an outcome of both the truncation itself and the choice for the 
initial binary fraction of 100%. Because of the lower mass limit 
the initial number of stars with m > 0.5 M increases for a given 
N sys , making the criterion on the halo more difficult to reproduce . 
In addition, since there is no truncation in the secondary mass 
distribution, a few VLMOs are part of a binary system and will 
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appeal" as single objects, either because their separation is larger 
than 400 AU or because the binary will be processed by dynam- 
ical evolution. For instance, in the case N sys = 20 and Rp\ =0.1 
pc, 2.5 VLMOs are identified in average at t = Myr. As a con- 
sequence the criteria on the VLMOs and on the halo remain dif- 
ficult to fulfil together with the other criteria. 
We ran additional simulations (model F) where we introduced a 
cut-off in the initial binary period distribution in a similar way 
as for model C (see section l67Zb . We find that when starting with 
a binary fraction of 40% and a cut-off at 50 AU many more sim- 
ulations could reproduce the six observational constraints with 
a success rate up to 10% for N sys = 20 and Rp\ = 0.1 pc. The 
average number of VLMOs in this initial configuration is 0.2 at 
t = Myr, which shows that the constraint on the number of 
VLMOs is easily satisfied, given that there are only close bi- 
naries (with separation < 50 AU) that very stable dynamically. 
It is interesting to note that the very few runs of model C that 
are also able to reproduce all the criteria correspond to the same 
initial conditions. Indeed the successful cases were obtained for 
A^ys = 20 and Rp\ = 0.05 or 0.1 pc and started with only one or 
two VLMOs initially. This seems to indicate that both the IMF 
and the initial binary population of tj Cha were not standard. 
When considering Murphy's constraint however, the success rate 
shrinks to 0.5% at best, if we require to have at least three stars 
in the mass range [0.08; 0.3] M and within a 10 pc radius. 
Despite a low success rate, this model is the only one that can re- 
produce all the observational constraints, including Murphy's re- 
sults. Note that the only successful runs are for two medium den- 
sity configurations: {N sys = 30; Rpi =0.1 pc} and {N sys = 30; 
fl P1 =0.05pc}. 

7. Summary and conclusion 

We have conducted a large set of pure N body simulations that 
aim to reproduce the peculiar properties of the tj Cha associa- 
tion, namely the lack of very low mass objects (m < 0.1 M ) 
and the absence of wide binaries (with a separation > 50 AU). 
We tested several models of various IMF and binary properties, 
and span the parameter space in density and virial ratio. The 
analysis was done using several procedures in order to compare 
efficiently the simulation results with the observational data and 
identify the best initial state. 

In order to test a universal picture for the IMF, we assumed a 
continuous log-normal single IMF with pi = 0.2 M and a = 
0.55. Starting with this IMF and a binary fraction of 100% (with 
either a random pairing, model A, or a flat mass ratio distribu- 
tion, model B), the analysis shows that ejecting all very low mass 
members without creating a halo of solar-type stars and keeping 
an inner core of 18 systems is not possible. Similarly to the case 
of a discontinuous single IMF, no simulation was able to match 
the observations. 

Reproducing all available observations of tj Cha by pure dynam- 
ical evolution from an universal single IMF and a stellar binary 
fraction of 100% is therefore very unlikely. 
We then tested a different set-up for the binary population, while 
preserving the shape of the IMF, our working hypothesis (model 
C). We assumed that wide binaries do not form initially and 
adopted a separation distribution truncated at large separation 
resulting in a lower binary fraction. As a result, the best initial 
state, starting with an initial binary fraction of 40% binaries and 
without any binary wider than 50 AU, yields a small success rate 
of 1 % (that drops to 0% if we requ ire those simulatio ns to have 
a halo of ejected low mass stars dMurphv et al.l |2010|) ). Since 
almost no considered initial state assuming a universal IMF sta- 



tistically matches the observational constraints, we started with a 
truncated IMF with no system below 0. 1 M e . However, this fails 
in reproducing the observations, unless starting with a singular 
binary population (no wide binary and a small binary fraction; 
model F). In this case, the best success rate is 10% and is ob- 
tained for initial parameters (N sys = 20 and Rp\ = 0.1; 0.05 pc) 
that are very similar to what is observed today in the cluster. 
This suggests that the dynamical evolution did not play a strong 
role in shaping the properties of 77 Cha and that most of them 
must be pristine, tj Cha may have started with an IMF deficient 
in VLMOs and with peculiar binary properties (namely a small 
binary fraction and an orbital period distribution truncated at 
small periods ) . Note t hat this conclusion is very different from 
iMoraux et al.l d2007bl) where the initial high density case was 
the preferred solution. This stresses the importance of the binary 
population in the overall dynamical evolution of the cluster. 
One can speculate onto the particular physical conditions that 
might have produced so few VLMOs together with preventing 
wide binaries from forming. 77 Cha may for instance originate 
from a hig hly magnetized cloud, preventing fragmentation of 
large scale dHennebelle et al.|[201 If) , forcing more mass into sin- 
gle fragments and not creating wide systems. Tighter binaries 
could then be produced later on, after the magnetic field has dif- 
fused out. 

Finally, in the low density case solution present ed above, it is 
very d ifficult to reproduce the recent results from Murphv et al. 
J2010I) . When considering this additional constraint, the success 
rate becomes very small (0.5% at best). Additional knowledge 
of the kinematics of this purported halo population might help 
refine the dynamical picture of tj Chamaeleontis. 
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